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Abstract. The spatial distributions of luminous and dark matter in massive early-type galaxies 
reflect the formation processes which shaped these systems. This article reviews the predictions of 
cosmological simulations for the dark and baryonic components of ETGs, and the observational 
constraints from lensing, hydrostatic X-ray gas athmospheres, and outer halo stellar dynamics. 

Keywords, (cosmology:) dark matter; galaxies: elliptical and lenticular, cD; galaxies: halos; 
galaxies: formation; galaxies: kinematics and dynamics; X-rays: galaxies; gravitational lensing 



1. Introduction 

A galaxy's dark matter mass is arguably the most important parameter characterizing 
its evolutionary state. Massive early-type galaxies (ETGs) with virial masses of up to 
~ 10 13 M0 are those that have evolved furthest in the galaxy building process, with little 
gas left to make new stars. In the centers of ETGs the luminous matter dominates, but 
already at a few half-light radii (R e ) of the stellar distribution, dark matter (DM) takes 
over. The spatial distributions of luminous and dark matter in the inner regions of ETGs 
are the outcome of the formation processes by which these systems were shaped, and 
are thus one of the observational constraints we have on these processes, next to scaling 
relations, stellar population properties, and orbital structure. 

The following sections review first the predictions of cosmological simulations for the 
structure of the dark halos and baryonic components in ETGs, and then what has been 
learnt about their dark matter halos from different techniques: lensing, hydrostatic hot 
X-ray gas athmospheres, and stellar dynamics of halo tracers. 



2. Dark matter halos from cosmological simulations 

2.1. Dark matter- only simulations 

Simulations of structure format ion in a Cold Dark Ma tter Universe have progressed 
greatly over several decades (see lDiemand fc Moorell2~011 , for a detailed review). For the 
idealized case of dark matter only, without baryons, resimulations of individual halos 
now reach O(10 6 ) particles per halo. The main results from these studies are: 

• Mass density profiles are approximately universal, i.e., independent of halo mass and 
cosmological parameters aside from scalings. The best current density profile model is 
a 3-parameter Einasto pr ofile ( Navarro et al. 19961: Bullock et al. 2001 ; Diemand et al. 
2004 iNavarro et al.|[2010h . 

• Halo concentrations are in the range c = r v { r /r s ~ 5-15 with scatter at given mass 
Alogc ~ 0.2. Early- forming halos and smaller halos have hi gher median concentrations 
( Navarro etaLlll996t iBullock et ailbOOll: iMaccio et al.|[2008h . 



• Halo shapes are strongly flattened, prolate-triaxial. The mean < c/a >— 0.5 ± 0.1. 
The elongation is larger at the center. It is approximately stable at fixed radius except in 
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major mergers, and is mostly due to anisotropy (median A c± 0.04 ) (jPubinski fe Carlbergj 
19911: Ijing fc Sutcll2002l lAllgood et aUbOOd: lHavashi et al.ll2007t) . 

• Bound substructures are abundant but not dominant, comprising ~ 10% of the 
halo mass. They have a steeply decl ining cumulative abundance distribu t ion with mass 

oc (V m ^/Vr 



max/ 'max.hnst I 



or OC 



M~ 2 (jKlypin et all Il999t iGao et all 12004 ISpringel et al 



20081 iDiemand et alJliooj ) 

• The inner density profiles are cuspy, increas ing inwards. The asymptotic inner slope 



is now known to be slightly shallower th an 
Dic mand et al ] l2004t iNavarro et alJlioioh . 



1 (Navarro et al. 1996; Moore et al. 1998 



2.2. Effects of baryons; alignment 

The structure of the simulated DM halos is modified by the gravitational interaction with 
the baryonic component, and thus indirectly by baryonic cooling, settling, star formation 
and feedback. The resulting modifications of the halos include: 

• Steepening of the ha lo density profile: adiabatic contraction through slow dissipative 
settli ng of gas into a disk (jBlumenthal et all 1986 ; Gnedin et al.ll2004 ; Sellwood fc McGaugh 
20051) . 

• Flattening of the cusp by feedback and shaki ng ( Gnedin fc Zhaol2002 ; Mashchenko et al 
20081 : iGovernato et alJlioiol : iMaccio et al.ll2012t ). 

• The halo b ecomes more spherical and evolves towards oblate at all radii but mostly in 
its inner parts feazantzidis et al. 2004 : Bailin et al. 2005 : Berentzen fc Shlosman 2006 : 



Abadi et al.l l2010) 



• Minor axes of the disk and inner halo (r < 0.1r V j r ) align when a disk forms. The mi- 
nor axis orientation of the outer halo remains unchanged. As a result, the orientations of 
the inner and outer halo are uncorrelated in simulations with baryons, whi l e well-aligned 



in N- body-only simulations ([Bailin et al 
20111) . 



20051 : iBailin fc Steinmetzl 120051 : iDeason et al 



The angular momenta of disk and inner halo also align, whereas the median mis 
alignment between inn er and outer halo angular momenta is ~ 25-45 deg (jBett et al 



20101 : lHahn et al.ll2010h . Because the halos rotate slowly, shape and angular momentum 
need not be tightly correlated. Stacking of simulated halos washes out the alignment 
signal for weak lensing. 

2.3. Mass structure of simulated early-type galaxies 
High-resolution models of ETGs in a cosmological setting are of particular i nterest for 



compa ring with the empirical results discussed in later sections of this review. IWu et al 
(|2012t ) analyzed 42 galaxies with stellar mass es of 2.7 x 10 10 -4 . 7 x 1 11 M from the 



cosmological hydrodynamical resimulations of lOser et al. (2010, |2012j) which include a 
model for cooling, supernova feedback, and star formation processes, but do not include 
AGN feedback. These systems contain an in-situ component made in efficient early star 
formation, and later grow by accretion through mostly minor mergers. 

Their stellar density distributions at z — are approximately described by Sersic 
profile s with realistic values of n and R e . As measure of their mass distributions. I Wu et al. 
(2012) determined circular velocity curves (CVCs) v c ee [GM(r)/r]^ 2 . The dark matter 



CVCs are essentially described by power laws outside radii greater than the softening 
radius (0.9/i _1 kpc for DM particles). They vary systematically in that the average 
slope of the dark matter CVC increases with increasing circular velocity, or mass, from 
approximately flat at low masses to slightly rising (~ 0.3) at high masses. When the stellar 
component is added in, the total CVC slope changes from slightly falling (~ —0.3) at low 
masses to flat at high masses (see Fig. la). The slopes of dark matter CVCs and total 
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CVCs correlate with n. Typical DM fractions in the simulated ETGs are ~ 20% at lR e 
and ~ 50% at 5 -R P .. They increase with R e in kpc, with n, or with mass within R e ; see 
Wu et al.l (|2012l ). These results are broadly consistent with the mass profiles determined 



from lensing, X-rays and dynamics, as described in the following. 



3. Dark matter halos of massive elliptical galaxies from lensing and 
X-rays 

3.1. Halo masses and shapes from weak lensing 

Halo masses for low redshift (z < .35) galaxies are constrained best by SPSS data (for 
other work with smaller sample s, see lBrainerd et al . 1996; iFischer et al. 2000; iHoekstra et al 
l2004l2005HGavazzi et al.ll2007t) . The weak galaxy-galaxy lensing signal for 350 '000 lenses 
(z < 0.35) from SDSS was compare d to a halo model bas e d on si mulations (dn/dM, NFW 
density profile, satellite model) bv iMandelbaum et al.l (|2006bl ). Satisfactory agreement 
was found; within this model an L* galaxy with M* = 6 x 10 10 M Q is hosted by a halo of 
mass 1.4 x 10 12 /i _1 M©. This value is similar for ETGs in high and low density regions. 
Results from the 2 ma g deeper RCS2 survey already suggested changes with redshift 
(jvan Uitert et al.ll201ll) ; see also the talk by Hudson in these proceedings. 

In principle, weak lensing is also the method of choice to measure the outer halo shapes: 
individual halos will have stronger shear signal on the projected major axis than on the 
projected minor axis. However, a measurement of this effect requires stacking of multiple 
galaxies. Thus the signal is preserved only if halos and galaxies are aligned. Simula- 
tions suggest t hat this is not t he ca s e and that the signal is isot r opized; see §2 . Lens - 
ing studies bv IHoekstra et al.l (|2004h : IMandelbaum et aD (|2006al ) : IParker et al.l (|2007t ): 
van Uitert et al.l ( |2012t ) have detected differing lensing anisotropy signals. These studies 



generally indicate triaxial halos but details are still unclear. 

3.2. Dark matter mass profiles from strong lensing 

Strong lensing (see, e.g., Treu 2010l ) gives an accurate measurement of the projected 
mass within the lens Einstein radius, Re- For lenses at z ~ 1, Re ~ several R e and is 
located where the DM mass dominates. In this case, the DM can be studied directly, at 
least in high-mass ETGs (a > 200 km/s) which dominate the samples. For most low- 
z lenses, Re < R e ] thus the strong lensing analysis measures the combined mass from 
stars and DM withi n Re- Mass profiles can be determined by combining mul tiple systems 
with different R E (|Rusin fc KochaneklEoOsT: iBolton et al.ll2008UGrilld|2012h or by mass 
recon struction methods applied to the source images ([Ferreras et al . 2005; ILeier et al 



20111 ) . The overal result from these studies is that on average DM starts to dominate at 



- 2-3i? e . 

By combining M(< Re) with a measured central velocity dispersion or resolved dis- 
persion profile, it is possible to constrain the average slope of the mass density pro- 
file, a. Results fro m the SLACS sample show that the density slopes are approximately 



isothermal (a = 2) (iTreu &; Koopmansll2004L Koopmans et al. 20061 ; Gavazzi et al.ll2007 



IBarnabe etal"1l2009f 



Auger et al.l l2010). This i s equivalent to a flat CVC in spirals and 



has been dubbed the "bulge-halo conspiracy" . IBarnabe et al. | (120091 1201 lh also measure 
halo axis ratios with axisymmetric models, obtaining values as flattened or rounder than 
the axis ratios of the stellar components. Adopting a Salpeter IMF, the DM fraction 
inside R E is found to increase fr om 25% fo r a = 200 km/s, M* = 10 11 M to 70% for 
a = 350 km/s, M* = 10 12 M Q (|Treull2010l) . Massive ETGs must ei ther be more dark 
matter dominate d in their cores, or have a heavier, non- universal IMF (jAuger et al.ll2010l; 
Treu et al.ll201ol) . 
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3.3. Hydrostatic mass distributions for X-ray bright ellipticals 

X-ray observations of the hot gas in massive ETGs have long b een an important and 
independent tracer of their DM halos (jBuote fc Humphrey! 120121) . The hot gas extends 
to large radii where the DM dominates, and the isotropic pressure makes the analysis 
much simpler than for stellar tracers. With Chandra and XMM-Newton, it is no longer 
a problem to measure temperature profiles. Thus the current main uncertainties are 
whether the gas is in approximate hydrostatic equilibrium, and how large are non-thermal 
contributions to the pressure which are not traced by the X-ray data, such as from 
relativistic particles, turbulent motions and magnetic fields. 
Some of the main conclusions are as follows: 

• In galaxies with X-ray luminosities Lx < 10 39 — 10 40 erg/s (depending on the optical 
luminosity, often lower mass systems) the X-ray flux from the hot gas is subdominant 
compared to that from X-ray binaries and s tars, so mass determination is unreliable 
(jRevnivtsev et alj|2008t iTrinchieri et al.ll2008l) . 

• A fr action of the more X-r ay luminous systems have strongly disturbed X-ray athmo- 
sphercs (iDiehl fc Statlerl | 2007h . so that the hydrostatic equilibrium assumption is doubt- 
ful. ICiotti fe Pellegrinil (|2004l ) quantified how inflows or outflows would disturb mass 
measurements. Reliable mass determinations are only possible in morphologically and 
dynamically relaxed systems (jBuote fc Humphrevll2012l ). 

• In X-ray bright and relaxed galaxies, the dark matter is clearly detected and the 
inferr ed potentials are near- i sothermal, i.e., the CVC is approximately flat to large 
radii ([Humphrey et all 120061 l2012t IChurazov et~aH |2008[ |201(A iNagino fc Matsushita 
20091). The inferred D M fraction at 2R e is ~ 40-80%. lDas et alj (f201oh using data from 



Churazov et al. ( 2010h found that the CVCs turn up at large radii, possibly due to the 



surrounding group DM halo, and that the typical circular velocity v c (R e ) correlates with 
the velocity dispersion of the sub-Mpc group environment. 

• Comparing determinations of the gravitational potential or CVC for the same galaxy 
from X-ray and optical observations gives a way to estimate the non-thermal contribu- 
tions to the total pressure. Values found are between ~ 10 % in quiescent galaxies like 
NGC 1399 to ~ 30% in more disturbed systems like M87 (|Churazov et al.ll2008l l2010t 



Das et al.l 120111 ). with large uncertainties because of possible systematic effects in both 
techniques. 



4. Dark matter halos from dynamics 

4.1. Kinematic halo tracers 

Traditional long-slit kinem atics reaches down to sur face brightnesses of Hv ^23.5 corre- 
sponding to R ~ 2R e (e.g. Kronawitter et al. 2000f) . To determine the DM mass profile 



and the halo orbit distribution, alternative data arc therefore needed which reach to 
larger radii and fainter surface brightnesses, such as 

• Planetary nebulae (PNe) velocities measured from their br i ght [OIII] emissio n line 
(e.g. iHui et al.|[l995l ; lArnaboldi et al.lll996t iMendez et all 1200 ll IPeng et alj|2004h. PNe 
are generally good tracers of the stellar light and kinematics in ETGs (|Coccato et al 
Useful kinematic information can typical! 



200' 



be derived to 
, reaching fxy r ~ 



- 5-8i? e 
27.5. 



and up to 



beyond 100 kpc in some cases (jDohertv et al.ll200S 

• Globular cluster (GC) velocities measured from their absorption line spect ra (e.g. 
ICote et alJ[200lt iHwang et all 120081: ISchuberth et alj|2010t IWoodlev et"al]|2010l) . Espe- 
cially bright, central ETGs contain large GC systems reaching to very large radii. The 
use of GCs as mass tracers is somewhat complicated by the fact that they do not trace 
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the galaxy light directly, requiring simultaneous determination of their (completeness- 
corrected) number density profile. 

• Stacked absorption line spectra from integral field units (IFU ) placed at large galac - 
tocentric radii. This method has be en employed using the Sauron ( Wei i mans et al. 2009T) 
and VIRUS-P (jMurphv et al.ll201ll ) IFU instruments and reaches down to fiy ~ 25.5. 

• Stacked absorption line spectra from off-ta rget pixels in slitlets originally placed i n 
the halo to measure globular cluster velocities ( Proctor et al. 2009t Foster et al. 2011 ). 
reaching [iy ~ 25. 

4.2. Outer halo kinematics from planetary nebulae 

While in early work with 4m telescopes velocities could be measured only for th e bright- 
est ~ 30-50 PNe in ETGs at Virgo distance (e.g. lArnaboldi et al.lll996l Il998l ). we now 
have ~ 20 ETGs with hundreds of measured PN velocities for tracing the velocity field 
out to large radii , typically 5-8R„. This is mostly thanks to slitless spectroscopy ef- 



forts at the VLT dMendez et al.ll200ll iTeodorescu et al.ll2005l l201ll : iMcNeil et al.H2oTol : 



McNcil-Movla n et al.ll2012D and to data f rom the special-purpose PN.S instrument (see 
Douglas et al.l 120071: ICoccato et al.ll2009t iNapolitano et al.ll2009[ l201ll). The kinemati c 
properties of ETG halos derived from these PN samples include (see Coccato et al.|[2009h : 

• The ratio of mean velocity v and dispersion a in the halo correlates with that within 
Re for the greater part of the sample, but some galaxies are more rotationally dominated 
in their halos. 

• The division into slow and fast rotators defined by the kinematics within ~ R e 
( Emsellem et al]|2007l ). is largely preserved in the outer halos. A few galaxies have more 
complex profiles of the angular momentum-related parameter. 

• Kinematic misalignments appear to be more frequent in the halos. 

• The radial profiles of rms velocity v 2 + a 2 fall into two groups, the major group char- 
acterized by a slow radial decrease, and a second group with "q uasi -Keplerian" , steeply 
falling profiles, first drawn attention to bv lMendez et al. ( 2001 ) and Romanowskv et al 
( 2003li . 



4.3. Kinematic substructure 

Cosmological galaxy assembly models predict that massive galaxies grow by accreting 
smaller companions. During and after minor merger events photometric and kinematic 
substructure would be visible in the outer halo of the host galaxy where dynamical time- 
scales are long. This is an interesting subject (and opportunity) by itself, but is also 
relevant for dynamical mass dete rminations at large radii. One examp le is the nearby 
BCG g alaxy NGC 1399, studied bv lMcNeil et~aT1 (|2010h with PNe and bv lSchuberth et al 



( 20101) with GCs. The phase-space plot for PNe around NGC 1399 in Fig. lb shows 
three components, including PNe from NGC 1399 around U1399 — 1400km/s, PNe from 
the lower-luminosity ETG NGC 1404 close-by in projection, and finally PNe from a 
component of "low-velocity outliers" with velocities ~ 750 ± 250km/s. 

The non-equilibrium low-velocity component is also seen in both the blue and red 
globular clusters around NGC 1399. The number density profile of red GCs approximately 
follows stellar light in this galaxy, whereas that of blue GCs does not. The decontaminated 
system of red GCs appears to be in approximate dyn amical equilibrium, but this seems 
not the case for the blue GCs (|Schuberth et alJ lioiO). This could be related to the fact 
that blue GCs are particularly prominent in smaller galaxies (jStrader et al. 2006) which 
would be preferentially accreted at late times. 

Indications for substructure have been foun d also in the halos of M87 and NGC 4649 
( Romanowskv et al.ll2012l ; Coccato et al. 2012 ). A fraction of large ellipticals have a sec- 
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ond galaxy within a 10-20 field. With discrete tracers (PNe or GCs) it is possible to 
eliminate kinematic substructures, which is a prerequisite for a reliable dynamical mass 
determination in such systems. 



4.4. Dark halo circular velocity curves and densities from stellar kinematics 

Using absorption line kinematics to l-2R e for luminous round ETGs, and nonparametric 
spherical distr i bution function models in combined luminous and dark halo potentials, 
iGerhard et al. ( 2001 ) found that the CVCs of ETGs are approximately fiat. A more recent 
analysis of a sample of ETGs in the Coma cluster with Mb =[-18. 8,-22.61 and kinematic 
data out to l-3i? e with axisymmetric Schwarzschild models by I Thomas et al. (2007) 
showed so mewhat more var ied CVC shapes. Massive ETGs tend to have fiat CVCs (e.g., 
NGC 4649. [Das et al. 2011 ). while lower-mass ETGs (see §4.5) are also consistent with 
slightly decreasing CVCs. These studies agree in their derived DM fractions, ~ 10 — 40% 
within R e , assuming constant M/L fo r the stellar component , which is also consistent 
with models of inner 2D kinematics bv lCappellari et all (|2006l ). 

Both studies also agreed in finding significantly higher central dark matter densities 
in ETGs than in spi ral galaxies of the same luminosity or mass. For the Coma ETGs, 
Thomas et al. (l2009h found a factor of 7x higher mean DM density within 2R e at the 



same stellar mass, and 13 x at the same luminosity. Baryonic contraction is not suffi- 
cient to explain this difference. One simple explanation is that the cores of ETG halos 



assembled earlier (a t redshifts 1-3) than spirals of same luminosity (jGerhard et al.ll2001 



Thomas et al. 1 2009) . Good agreement of t he dynamical l umino us p lus dark matter model s 



was found with SLACS lensing models of Auger et al.l ( 20091 ): see lThomas et al. ( 2011 ) 



4.5. Mass Distributions of Quasi- Keplerian Ellipticals 



Mendez et al.l (|200ll ) first showed that the steeply decreasing outer velocity dispersion 
profile of the intermediate-luminosity elliptical (ILE) galaxy NGC 4697 could be matched 
well with an isotropic Hernquist model. They thus concluded that no evidence for dark 
matter had been found out to 3R e , but that dark matt er could be present if the velocity 



distribution is anisotropic. Romanowskv et all ( 2003 ) drew attention to three further 
ILEs with steeply decreasing a p (R), NGC 821, 3379, and 4494, suggesting the presence 
of little if a ny dark matter in their halos. However, the well-known mass-anisotropy 
degeneracy (jBinney fc Mamonl Il98i ) is much stronger in galaxies such as these ILEs, 
with de Vaucouleurs type luminosity profiles and steeply decreasing er p -profiles, than 
in galaxies with eit her more shallo w luminosity profiles or radially constant projected 
velocity dispersions (lGerhardlll993h . 



With NMAGIC particle models (see lde Lorenzi et al.ll2007l ) ba sed on a variety of kine 



matic data, incl uding several hundred PN line-of-sight velociti es. Ide Lorenzi et al.l (12008 



for NGC 4697). Ide Lorenzi et al.l (|2009l for NGC 3379) and iMorganti et al.l (|2012l for 



NGC 4494) determined the range of quasi-isothermal halo potentials consistent with the 
data for these galaxies. Fig. lc shows the 70% confidence boundaries from the PN like- 
lihood for NGC 3379 and 4697, and from a more comprehensive analysis for NGC 4494. 
While some spread in the allowed CVCs and anisotropies remains, as expected, dark 
matter is required in all 3 galaxies. Fig. lc indic ates that the mos t favoured CVCs may 
be slightly falling, as in the lower-mass systems in Wu et all ( 2012t ). Relative to the outer 
CVCs, the baryonic centers of these galaxies are quite centrally concentrated, as could 
be expected in gas-rich mergers. 
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Figure 1. (a) Top: Circular velocity curves from stars (magenta), dark matter (black-blue, 
both left pan el) and total mass distribution (right), for 42 model galaxies from the cosmological 
simulation of lOser et al.l 1120121), in b ins of increasing total circular velocity at 5R e , as given 
on the right panels. From IWu et alj (|2012f ). (b) Bottom left: Phase -space distribution of PNe 
in the nearby BCG galaxy NGC 1399 (black), its neighb our galaxy NGC 1404 (blue), and 
the low- velocity component (red). From Mc Neil et alj (|2010l ). (c) Bottom right: 70% confidence 
ranges from dynamical modeling for total circular velocity normalized to value at R e , DM 
fracti on, and anisotropy par ameter, for the three quasi-Keplerian ETGs NGC 3379, 4494, 4697. 
From lMorganti et al.1 (|2012f ). 
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4.6. Dark Matter Fraction with Mass 

For radially constant mass-to-light ratio T, maximizing the fraction of mass following 
the luminosity profile results in a lower limit for the DM fraction, /DM,min- Dynamical 
modeling (§4.4, 4.5) gives / D M,mm(-R < IRe) ^ 0.1-0.4 and / D M,min(i? < 5i? e ) ~ 0.3-0.8, 
with a tendency of larger valu es for more massive E TGs. The corresp o nding T max also 
increase with galaxy mass fsee lThomas et al1l201l[ ). iNapolitano et al.l ( 20051 ) estimated 



gradients of T for a sample of ETGs with available dynamical mass modeling and found 
a correlation between luminosity and T-gradient such that the more luminous galaxies 
were more dark matter dominated. Comparing with DM halo predictions, the underlying 
reason appeared to be that the stellar bo dies of brighter galaxies enclose a larger fraction 



of their DM halos. iDeason et al.l (|2012f) used power-law models to estimate dynamical 
masses at 5i? e for a sample of ETGs with extended PN and GC tracer kinematics. From 
their models they also found an increase of DM fraction inside 5R e with galaxy mass, 
independent of whether a Chabrier or Salpeter IMF was assumed. The strong lensing 
results (§3.2) show a similar trend for the DM fraction within Re- 

Because dynamical T max increase with galaxy mass faster than stellar population T 
for given IMF, either some of the dark matter in massive ETGs must follow the ligh t 



2011). 



distribution, or the IMF must vary systematically with galaxy m ass (I Thomas et~al 
The latter interpretation is favoured by the dynamical analysis of lCappellari et al.l (|2012l ) 
and by new stellar populati on modeling of ETG spectra with NIR spectral features 



spectra 

sensitive to low mass stars ( Conrov fc van Dokkumf[20rl . If the IMF varies between 



galaxies, it is natural in a hierarchical assembly picture to expect that it would also 
vary with radius within a galaxy. Then reliable mass decomposition in ETGs will require 
radially resolved spectral measurements of the IMF. 



5. Conclusions 

• Pure N-body simulations in ACDM cosmology predict nearly universal, prolate- 
triaxial DM halos. Including baryons and galaxy formation modifies the DM halo shapes, 
alignment, and inner density profiles. This provides one way to constrain the formation 
processes of the baryonic components. 

• Weak lensing measurements are consistent with the predicted density structure of 
halos on large scales, and determine the relation between halo mass and galaxy luminosity. 

• Mass determinations from strong lensing, X-ray emitting gas, stellar kinematics and 
PN kinematics out to several R e find flat, close to isothermal CVCs for massive ETGs, 
while for lower mass, "quasi-Keplerian" ETGs the CVCs may be slightly falling. 

• DM fractions at 5i? e inferred from dynamical modeling are in the range ~ 30-80%. 
Inner DM densities in ellipticals are ~ 10 times higher than in spiral galaxies of the same 
stellar mass. 

• A suite of recent high-resolution cosmological simulations of massive galaxies grown 
by minor and major mergers show systematic variations of the CVCs and DM fractions 
inside 5i? e with mass and fitted Scrsic index, similar to those observed in ETGs. 
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